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Abstract: Using numerical ray tracing, the paper studies how the average distance modu- 
lus in an inhomogeneous universe differs from its homogeneous counterpart. The averaging 
is over all directions from a fixed observer not over all possible observers (cosmic), thus 
is more directly applicable to our observations. In contrast to previous studies, the aver- 
aging is exact, non-perturbative, and includes all non-linear effects. The inhomogeneous 



universes are represented by Swiss-cheese models containing random and simple cubic lat- 



o 

tices of mass-compensated voids. The Earth observer is in the homogeneous cheese which 
f*T} [ has an Einstein - de Sitter metric. For the first time, the averaging is widened to include 

the supernovas inside the voids by assuming the probability for supernova emission from 
any comoving volume is proportional to the rest mass in it. Voids aligned along a certain 
direction give rise to a distance modulus correction which increases with redshift and is 



o 



caused by cumulative gravitational lensing. That correction is present even for small voids o 
and depends on their density contrast, not on their radius. Averaging over all directions 
0^ . destroys the cumulative lensing correction even in a non-randomized simple cubic lattice 

of voids. At low redshifts, the average distance modulus correction does not vanish due to 
the peculiar velocities, despite the photon flux conservation argument. A formula for the 
maximal possible average correction as a function of redshift is derived and shown to be in 
excellent agreement with the numerical results. The formula applies to voids of any size 
that: (1) have approximately constant densities in their interior and walls; and (2) are not 
in a deep nonlinear regime. The average correction calculated in random and simple cubic 
void lattices is severely damped below the predicted maximal one after a single void diam- 
eter. That is traced to cancellations between the corrections from the fronts and backs of 
different voids. The results obtained allow one to readily predict the redshift above which 
the direction-averaged fluctuation in the Hubble diagram falls below a required precision 
and suggest a method to extract the background Hubble constant from low redshift data 
without the need to correct for peculiar velocities. 
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1. Introduction 

It is fair to say that the standard cosmological ACDM model is facing a phenomenologi- 
cal crisis. The dark matter carrier has been evading direct detection for decades and the 
origin of dark energy remains a theoretical puzzle. The most natural candidate is the 
vacuum state energy, but the flat-space Quantum Field Theory is incapable of calculating 
it, producing an estimate that is 10 120 times as large as the value suggested by the su- 
pernova observations. Presumably, that would be rectified in a fully quantized theory of 
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gravitation which unfortunately does not yet exist. Remaining within the standard Gen- 
eral Relativity, there are attempts to explain dark energy as an apparent quantity arising 
from the averaging procedure that maps the real lumpy universe to the idealized homo- 
geneous Friedman-Lemaitre-Robertson- Walker (FLRW) metric. Such approaches separate 
into three distinct classes. First, one could average the Einstein equations over spatial 
hypersurfaces and hope the analogue of the Friedman equation contains a significant ef- 
fective A term. Such an idea is conceptually problematic since inhomogeneous universes 
do not naturally pick up preferred spatial slices to average over, unlike the homogeneous 
FLRW models. It is shown in Q that the result of such averaging depends on the arbitrary 
choice of slicing and has no coordinate independent meaning - even Minkowsky spacetime 
could produce an apparent acceleration. A cousin to the first approach is averaging over 
the null hypersurface of the past light cone of the observer, g, ||. It is not proven that 
such averages have a physical meaning and are not simply coordinate quantities. In the 
second approach, the second order perturbations to the Einstein equations are calculated 
and treated as an effective energy momentum tensor term. Unfortunately, as pointed out 
in the second order perturbation term is not gauge independent and therefore is not 
suitable to play the role of an energy momentum tensor in the Einstein equation. 

The present paper belongs to the third approach which is completely coordinate inde- 
pendent since it calculates relations only between physically observable quantities such as 
redshift and luminosity distance (or its logarithmic measure - the distance modulus). Pa- 
pers of this type, for example demonstrate that the distance modulus - redshift Hubble 
diagram of type la supernovas can be reproduced without any dark energy for an observer 
occupying the center of a giant (~ Gpc) underdense spherical bubble. That is possible be- 
cause the central observer measures a local Hubble parameter that is larger than the global 
one. The corresponding shift on the Hubble diagram between the true background and the 
one assumed by the observer allows for fitting the supernova data. Bubbles of such a large 
scale have significant peculiar velocities away from their center. That creates a fine tuning 
problem - the observer has to be in the vicinity of the bubble center || to observe a small 
CMB dipole equal to the measured v/c ~ 0.002 (Local Group velocity 620 km/s). This 
violates the Copernicus principle that we do not occupy a special position in the universe. 

Other papers of the same class studied the collective effect of a configuration of 
many bubbles/voids of smaller size, thus avoiding big peculiar velocities and significant 
anisotropies in CMB. That scenario can be conveniently simulated in the "Swiss-cheese" 
toy model j|, 0, ^] which has the virtue of being analytically solvable. It is constructed 
by removing spherical regions from a homogeneous FLRW background (the "cheese") and 
replacing them with inhomogeneous density distributions with the same gravitating mass 
(mass-compensated voids). It was used in |], |3| for an observer looking along the diameters 
of a string of aligned voids of radius 350 Mpc. A positive cumulative change in the dis- 
tance modulus was found which increased with redshift with respect to the homogeneous 
background. Although the size of the correction was significant at high redshift, it was 
not large enough to fit the supernova Hubble diagram and it substituted only partially for 
dark energy. A smaller correction due to aligned voids with a more mundane radius of 
23h~ l Mpc was calculated in [1C]. The result obtained in || was criticized in [11] where it 



was demonstrated to stem from the cumulative weak lensing defocusing of the rays passing 
diametrically through the aligned voids. The same paper showed that the correction to 
the distance modulus vanishes when averaged over randomized impact parameters of the 
rays entering the voids. This conclusion is in accord with an old argument |]l^] based on 
the gravitational lensing conservation of the total photon flux. The implicit geometrical 
assumptions in \12\ have been challenged in more recent papers ]i~3| , |Ti| ] (and references 
therein) and the present work will demonstrate they are violated close to the observer due 
to peculiar velocity modifications of the redshift surfaces. 

The studies mentioned above are deficient in several areas. Averaging over directions 
was not performed in Q. The calculations in were done using weak lensing theory 
neglecting the time dependence of the void density and possible nonlinear effects which 
become significant at low redshifts. The averaging assumed that the ray impact parameters 
had a uniform probability distribution over the cross-section of a void. That is increasingly 
true at high redshifts but does not hold in the local neighborhood where the rays have to 
converge on the observer. More importantly, [11| averaged only over supernovas residing 
in the homogeneous cheese. The same was done in a study that obtained an analytical 
estimate of the corrections to the luminosity distance, [jig ]. Such a bias is unjustified, 
first because supernovas in the voids produce much bigger corrections than the ones in the 
cheese ||, and second, the observed supernovas occur more frequently in denser regions 
like the void walls than in the cheese. 

The goal of the present paper is to address once again the question whether the 
direction-averaged distance modulus correction from a configuration of voids really van- 
ishes but this time in the context of a calculation that: (1) is exact, non-perturbative, an 
includes all possible non-linear effects; (2) is armed with a physically sensible averaging 
procedure free of ad hoc assumptions about impact parameters [11] or random cancellations 
of corrections fl2|| ; and (3) does not neglect the supernovas inside voids. 

The final outcome of the calculations cannot be guessed on the grounds of the previous 
studies because it may depend on non-linear effects, the choice of averaging procedure, and 
the supernova selection, none of which was taken into account before. The calculations 
are performed in Swiss-cheese models with mass-compensated voids of two radii having 
observational support: 30 and 300 Mpc. The cheese is spatially-flat matter-only Einstein 
- de Sitter (EdS). The observer is placed in the cheese since at present there is no indi- 
cation that we occupy a void. The first Swiss-cheese model is set up in section 2. The 
observer shoots past-directed light rays in all directions. The light propagation geodesic 
equations and the luminosity distance tracing along each ray are discussed in section 3. 
The physically sensible averaging procedure in this paper assumes that the probability 
for supernova emission from a comoving volume is proportional to the rest mass in it. 
Averaging the distance modulus for a single void is discussed in section 4 and it carries 
the essential characteristics of the procedure for many voids. The main outcome of that 
section is a simple formula to estimate the maximal average corrections to the distance 
modulus. Section 5 studies the cumulative correction due to gravitational lensing along a 
string of aligned voids. Numerical results from averaging within random and simple cubic 
void lattices are presented in section 6. The effect on the distance modulus produced by 



large voids of radius 300 Mpc, which leave a measurable imprint on CMB, is evaluated in 
Section 7. The summary and conclusions section discusses the answer to the main question 
addressed by the present study and the practical importance of the obtained low-redshift 
results for future surveys. 

The speed of light in this paper is c = 1, and times and distances are measured in 
megaparsecs (Mpc), occasionally giving time in megayears (Myr) for convenience. The 
gravitational constant G is kept explicit in all equations so that the reader can easily 
substitute with a favorite value. The usual geometrized units, G = 1, were used in the 
numerical calculations. 



2. Swiss-cheese model with voids of radius r v = 30 Mpc 

2.1 Einstein-de Sitter as homogeneous cheese 

Every reasonable model of the universe is matter dominated in the past while the structure 
is being formed. Adding a cosmological constant has little impact [^] on the development 
up until recent times when it starts dominating the matter. That motivates choosing the 
spatially-flat matter-only Einstein-de Sitter (EdS) model for the homogeneous regions of 
the Swiss-cheese. The matter includes both visible and dark varieties. This is a convenient 
playground to explore the question whether voids in the matter distribution can mimic 
the effect of dark energy. The Hubble constant of the model is Hq = 70 km/s/Mpc. 
The big bang time is conventionally chosen as tbb = 0, the current age of the model is 
to = 2/(3ifo) = 2857 Mpc (9312 Myr) and the scale factor as a function of cosmic time is 
a(t) = (t/t ) 2/3 . The matter density today is p(t ) = 3H$/(8nG) and as a function of time 
is: 

p{t) = Wf = ^gT 2 ■ ( } 

2.2 Lemaitre model of a spherical void 

The voids in the Swiss-cheese are modeled with the Lemaitre metric |], [7j also known as 
" Lemaitre- Tolman-Bondi solution". It describes a spherically symmetric spacetime filled 
with an irrotational pressureless ideal fluid (matter, dust). The matter particles are in a free 
fall under their own gravity tracing geodesies. The zero rotation of the geodesic congruence 
means the geodesies are hypersurface-orthogonal and the corresponding family of spatial 
hypersurfaces define a convenient foliation and coordinate system on the spacetime. The 
resulting coordinates = (t, r, 9, </>) are matter-comoving and synchronous. The matter 
energy-momentum tensor in these coordinates is diagonal: T^ u = diag(p(r,t), 0, 0, 0). The 
metric is given by: 

E' 2 

ds 2 = -dt 2 + — — — -cir 2 + R(r, tfidO 2 + sin 2 9 # 2 ) , (2.2) 

1 + 2E[r) 



where prime denotes a derivative with respect to the radial coordinate r. The arbitrary 
integration function E(r) results from integrating the Gt r = Einstein equation. The areal 
radius R(r, t) determines the area of a sphere of radius r and E(r) determines the local 



3-curvature of the spatial slices. The time coordinate t measures the proper time of the 
comoving matter. Integrating once the G rr = Einstein equation leads to the following 
evolution equation: 

■ n . , 2GM(r) , . 

where dot denotes a derivative with respect to t and M(r) is another arbitrary integration 
function. The above equation alludes to conservation of kinetic plus gravitational energy 
in Newtonian mechanics where the "energy function" E(r) plays the role of total energy 
per unit mass. The rest of the Einstein equations connect the "mass function" M(r) to 
the comoving matter density p: 

M'(r) C r 

" (r - 4> = • M M " I *"*'- • f t)if ' < 2 ' 4 ' 

The integral can be evaluated for any t and will produce the same M(r). Note, M(r) differs 
from the comoving rest mass which has an extra factor of (1 + E(r))~ l l 2 in the integral. 
This is an example of a "relativistic mass defect", M, M. The parametric solution of ( |2.3| ) 
for the E > case of interest is: 

GM(r) GM(r) 
R{r > t] = ~2eU (C0Sh7? " 1} ' * " tB{r) = (2^(r))3/2 (siDh7/ " V) ' (2 ' 5) 

where the "big bang time" is another arbitrary function and rj > is a parameter. 

It is natural for inhomogeneous models to have different places with big bang happening 
at different times. 

The metric fl2,2|) contains as a subclass all FLRW models which are obtained by setting 
R = a(t) r and E = —kr 2 , k = const. In particular, the cheese regions of the present model 
are spatially- flat and have E = 0. The voids have to be matched to the homogeneous EdS 
without tearing of the metric. The appropriate junction conditions 0] are: 

E{r v ) = 0, t B (r v ) = 0, M(r v ) = M EdS (r v ) , (2.6) 

where r v is the radius at which the void merges with a homogeneous region. The first two 
conditions express an obvious continuation to the EdS values. The third condition means 
the void is "mass-compensated": the gravitating mass inside it matches the mass in EdS 
within the same radius r v . The conditions above guarantee that the metric outside radius 
r v remains EdS exactly - it does not "feel" the presence of the void and the metrics inside 
different voids do not interfere with each other. Of course, voids in the real world are not 
restricted within some radius, they interfere and even merge, and their walls accrete mass 
from outside and regions where the metric stays exactly homogeneous do not exist. The 
inconvenience is that a real world simulation requires computing of the global metric (or 
its Newtonian approximation) encompassing the whole spacetime. The mass-compensated 
voids are a toy-model which is easier to compute because the metrics inside different voids 
are exact copies of each other. 

The Lemaitre metric is specified by choosing three functions, say R(r,ti), E(r), and 
M(r), where is some initial time and the last two functions determine R(r,t) via (|2.3|). 



The radial coordinate r is just a label for the spherical matter shells and can be chosen 
at will (gauge freedom). That means the system is specified by only two functions plus 
a gauge choice for r, effected by choosing R(r,ti). Various possible selections of the two 
functions are discussed in [17|. 



The r gauge in the present paper is fixed by choosing 

R(r,t ) = r. (2.7) 



It was demonstrated in [18] that the Lemaitre model can be written from synchronous to 



perturbed FLRW coordinates (known as "Newtonian gauge"): 

ds 2 = -(1 + 2ip)dt 2 + a(t) 2 (l - 2ip){dx 2 + dy 2 + dz 2 ) , (2.8) 
as long as the peculiar velocities remain small (which is evident later). The void model 



needs initial conditions in the synchronous coordinates (2.2) but the astronomical data 
uses the Newtonian coordinates ( |2.8| ). This dichotomy is resolved by the choice in Q2.7| ) 
which makes the two coordinate systems coincide, approximately, at time to [jig]. That 
allows for using Newtonian notions in synchronous coordinates around time to an d feeding 
astronomical data of distances, densities and velocities directly to the model. Other models 
in the literature fix the synchronous r gauge at some initial time much earlier than to- By 
the time to, their synchronous coordinates evolved significantly away from the Newtonian 
ones which prevents direct Newtonian interpretations and may lead to illusory coordinate 
effects. For example, light rays that look straight in Newtonian gauge often appear "curved" 
in such synchronous coordinates. 

The present paper uses p(r,to) and ts(r) as the two functions specifying the void 
metric. There is not a standardized definition of what constitutes a void in astronomy. 
The so called supervoids are defined as regions free of rich clusters and in the Milky 
Way proximity have an average radius of about r v ~ 46/i _1 Mpc ~ 65 Mpc (for Hq = 
70 km/s/Mpc) |I~9[ |. Computer void-finding algorithms usually define a void nucleus to 
be completely free of galaxies leading to breaking down the supervoids into smaller ones 
of r v ~ 20 Mpc PQ|. Different algorithms often find voids of different sizes, shapes and 



numbers in the same region of space [21|. A recent paper examined the SDSS 5th release 



data |22| using an algorithm that finds voids closest to visual inspection and found an 



average r v ~ 36 Mpc [22]. It is obvious there is not an universal agreement upon the 
average void size and the current paper adopts an intermediate value of r v = 30 Mpc. 
The conclusions reached at the end scale with r v as long as the voids are not significantly 
nonlinear. The current galaxy density inside voids is very low, below 10% of the average. 
Density reconstruction of our neighborhood using the peculiar velocity field shows that the 
total density (dark and visible matter) inside voids drops below 0.40 p(io) [ |23| ] . The bias 
between visible and dark matter is still an open question but the present paper will assume 
that light is a good tracer for matter and set the total matter density inside the void to 
Pin(to) = 0.10 p(to). Observationally, most of the visible matter is gathered in the void 
walls and it will be assumed the same applies to the dark matter distribution as well. The 



wall thickness is set as ~ 5 Mpc. The so designed matter density of the void model is 

29)] 



p(r,t ) = p(t ) 



Ax + A 2 tanh[2(r - 25)] - A 3 tanh[3(r 
1 



2.9) 



,r < 30 
,r > 30 

and is shown on FigJTJ. The values of the constants to six significant figures are (A\, A2, A3) 
= (0.548006, 1.25446, 0.806455). They were determined by requiring density continuity 
and satisfaction of the mass-compensation condition in (2.6). The mass function is obtained 
by evaluating the integral in (2.4) at time to- 

The Lemaitre metric is fixed completely by specifying one more function. This could 



be the peculiar velocity at time to, [17], but it is not clear what a typical velocity profile 



looks like. The bang time t B (r) turns out to be a better candidate. After specifying it, 



one needs to solve for E(r) the two equations in (2^5) evaluated at time to. Eliminating 77 
gives 

arccosh{l + Aqx) = x 3/2 (t - t B (r)) , (2.10) 



where Ao(r) 



((i + A)z) 2 -i) 1/2 

= R(r,t )/{GM) 1 /' 3 and x(r) 



(2.10) to have a solution x(r) at a given r is derived in JTj 

t B (r) > t crit (r) = t 



2E/(GM) 2 / 3 . The condition for equation 



(2.11) 



Equation (plpp can be solved numerically for x(r) which in turn gives E(r). All the models 
with the given density ( [2.9D are enumerated by the possible functions t B (r) > t cr u(r). The 
present model sets 

t B (r) = 0, (2.12) 

which satisfies the junction condition (|2.6| ) and the solvability condition (|2.11| ). 

Choosing dt B /dr = excites only the p < r , r )/p(r ) 
growing density perturbation mode |24], p5| 
which is favored by linear perturbation the- 
ory in cosmology and increases like 5 = 8p/p 
oc t 2/3 oc a(t) in EdS |2§. The decaying 
mode 5 oc H{t) = a/a [27] decreases like 
5 oc t _1 in EdS. It is neglected in cosmol- 
ogy assuming it had enough time to decay 
to insignificant levels. Choosing dt B /dr 7^ 
excites a mixture of growing and decaying 
modes. A common example are Lemaitre 
models that set the peculiar velocity at ini- 
tial time ti to zero which in effect creates a sum of growing and decaying modes with 
amplitudes 5 grow (ti)/5 decay (ti) = 3/2 |2§. 

The model selected by ( p,12| ) is a pure growing mode. Knowing E(r), the metric 
function R(r, t) can be obtained by either using the analytic solution (|2.5| ) or by numerically 
integrating fl2.3|). The second method is advantageous in Mathematica since it produces 
an interpolation function which is fast to call and whose numeric derivatives are easy to 
calculate. 




r, Mpc 



Figure 1: Current void density p(r,to) divided 
by the homogeneous density p{to). 



2.3 Void model properties: density, shell crossing, last scattering, peculiar ve- 
locity 

The comoving matter density of the void is shown on Fig.|2] and behaves as expected. The 
ratio of the void wall density to the homogeneous density is increasing at all times. The wall 
eventually collapses, reaching infinite density at shell crossing which happens at t « 4500 
Mpc (14680 Myr). 



p(r,t) / p(t) 




P (r, t m ) I p (r LS ) 
l.OOlf 

1 













r, Mpc 







20 25 30 35 



r, Mpc 



Figure 2: Density ratio p(r,t)/p(t) at times 
(500, 1000, 2000, t , 3500, 3900) Mpc (from 
bottom to top). 



Figure 3: Density ratio p/p at time of last 
scattering. 



The time of last scattering in this model, defined as the time at which the EdS scale 
factor is a = 1/1100, is tis = 0.0783 Mpc (0.255 Myr). The density ratio at last scatter- 
ing is shown on Fig.||. The density contrast of the perturbation is of the expected order: 
5 = p/p — 1 ~ 10 -3 . This 6 applies to the dominating dark matter; the baryonic mat- 
ter perturbations, suggested by the CMB data, are still 5bar ~ 10 -5 at the time of last 
scattering since the baryons had just decoupled from the photons. 

The peculiar velocity is defined as the 
matter velocity with respect to the Newto- 
nian coordinates (|2.8|). This concept does 
not apply to synchronous coordinates at all 
but can be calculated from quantities in them, 
@: 

a(t) 



v(r, t) 



0.0020 



0.0015 



0.0010 



0.0005 



v(r,t) = R(r,t) - R(r,t) 



a(t) 



v < 1, 
(2.13) 




t- r, Mpc 



Figure 4: Peculiar velocity profiles at times 
t=(t LS , 10, 100, 1000, t , 4500) Mpc (from bot- 
tom to top). 



Fig.|] shows a few velocity profiles at different 
times. The peculiar velocity is everywhere 
positive since in Newtonian coordinates ( |2.8| ) 
the matter underdensity inside the void acts 

like a repelling gravitational source which accelerates the matter at bigger radius outwards. 
Even if the initial velocity was negative (towards the center), it will be reversed at later 
times outwards. For an underdense void interior, a negative initial velocity there is possible 
only if the model contains a decaying perturbation mode. The current model has a pure 
growing mode and the peculiar velocity is always positive. 
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The velocity peak in the profiles can be estimated with standard linear perturbation 
theory. The linearized mass continuity equation in Newtonian coordinates gives, p8|j . 
V • v = —ad5/dt = —a5 din 5 /din a. If the void interior is homogeneous, the velocity 
there is linear (see Fig.[|): v(r,t) = n(t)r, where (r, t) are Newtonian coordinates fl2.£| ) 
comoving with the averaged matter, not the synchronous ones. For the growing mode, 
8 (x am. the EdS. Plugging back in the continuity equation obtains 3k = — a 5 which leads 

t0 11 

v = — ar5 = — Hard. (2.14) 
3 3 y 1 

Fig.|] shows that the velocity peaks at synchronous r ~ 25 Mpc. Taking that value for the 

Newtonian r and noting that today 5 (to) = —0.9 in the void interior one gets v max (to) ~ 

1.75 x 10~ 3 when the actual value from Fig|4] is v max (to) ~ 2.26 x 10~ 3 . The agreement is 

not bad if one notes that the perturbation is already in a nonlinear regime at time to- 

The shape of the velocity profiles is explained by the velocity asymptotic behavior 

close to the big bang when t — tg(r) — > and rj — >■ 0. One can expand the second equation 

in (2J3) and then invert the series to obtain 

' = 6l,3B -To B3 + ^k B * + ° (b7 > ' B( - ( > s ^> ( w) " 3 ■ (2 ' 15) 

Substituting that in the first equation gives 

R f r , t) = ^r\ ^B 2 + -t^B A - —B 6 + 0(B 8 )) . (2.16) 
K ' 1 2E 12V3 2 2 / 3 10 1400 y 'J v ' 

Inserting the above equation into the peculiar velocity definition fl2.13| ) and using a/a = 
2/(3t) for EdS obtains 

t-t B y/ 3 ((A/^Hb , (3/4) 1 /3 ( i + t B / t ) g^ + ts/t)^ , „ f A 



(2.17) 

The first term in the above expansion has a time dependence (t — tg) 1 ^ tg/t ~ t -4 / 3 for 
small tg <C t. One can identify it with the decaying density mode, |2q| . As mentioned 
before, that mode vanishes for models with tg = const (=0 by (|2,6j )). In this case, for 
times not too far in the future: t <C GM/(2E) 3 ^ 2 (ensuring B < 1) , the second term 
in the above expansion dominates and the velocity profile evolves in time as oc t 1 / 3 . This 
corresponds to the velocity of the growing density mode, pq] . The constant shape of the 
velocity profile is set by E(r)/(GM(r)) 1 / 3 which is typically triangular since that quantity 
is zero at r = and r = r v and positive in between for E > models. 

3. Ray tracing 
3.1 Redshift 



Light rays are null geodesies. The 4-momentum of a photon along such a geodesic inside a 
void is = ko (dt/dX, dr/dX, d9/dX, dcp/dX), where ko is a constant depending on the 



choice of the affine parameter A. The redshift z of a photon emitted by a supernova and 
observed on Earth satisfies 

-. . Kbs u em 

= (U ■ K) em = (dt/dX)em (3 x) 

Km ^obs (U ■ K) obs (dt/d\) obs ' 



where is the 4- velocity of the emitter and the observer in the comoving coordinates (£ 
and dot denotes scalar product. Both the supernova and the Earth observer are assumed 
comoving with the matter therefore = (1, 0, 0, 0) which leads to the last equality. The 
affine parameter is uniquely defined up to scaling and shifting (affine transformations). 
One can use those to choose a convenient parameter that satisfies 

Kbs = 0, (dt/d\) obs = 1 , (3.2) 



where obs means evaluated at the Earth observer. Substituting that in (fO|) gives the 
connection between time and redshift along a light ray: 

3.2 Local and global spatial coordinates for ray tracing 

The affine parameter A is chosen to increase with the time t along the ray. The ray's 
geodesic equations are best integrated by following the ray back in time (decreasing A) 
from the Earth observer (A = 0) to a supernova (A < 0). Integration forward in time would 
require knowing the precise ray direction at the supernova required to hit the Earth and 
that is hard and impractical to calculate. Note that tracing a ray by decreasing A does not 
flip the sign of any derivative with respect to A! 

The fastest method of ray integration inside a void is to use local polar coordinates 
(r, 9) in the two-dimensional plane containing the ray and the void center, Fig.||. Due 
to the symmetric metric (|2.2[) , a light ray always remains in that plane which can be 
taken to be <f> = const by properly orienting the spatial coordinate axes. This reduces the 
number of independent variables by one and simplifies all equations. From void to void, 
the rays are traced in global spatial coordinates comoving with the homogeneous matter. 
In the homogeneous cheese, these coordinates are the usual FLRW Cartesian 3-vectors 
x = (X, Y, Z). Inside a void (for r < r v ), the global coordinates are Cartesian-like and are 
connected to the local (r, 9) by the familiar relation: 

x = x c + r cos 9 u + r sin 9 v , (3-4) 

where x c contains the global spatial coordinates of the void center C and the unit 3-vectors 
{u, v} form a basis in the two-dimensional plane of the ray, Fig.||. The above equation 
does not carry the usual physical sense of conversion from polar to Cartesian coordinates - 
the spatial metric inside a void is not really Cartesian if written in the global coordinates 



(|3.4| ). Despite that, the coordinate transformation ( |3.4[) is mathematically permissible and 
one is free to use it. 

The u vector is chosen to point from the void center to the point F at which the ray 
leaves the void (enters the void when propagating back in time). The angle 9 is measured 
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Figure 5: Local coordinates in the ray plane: C - void center, F - point of ray exit, {u,v} - basis 
unit vectors. 

from the u-axis. The v vector is chosen orthogonal to u pointing in the half-plane occupied 
by the ray. The global components of {u, v}, needed to calculate the components of x, are 
determined geometrically from the global coordinates of the void center, x c , the point F, 
and the ray direction at F. 

A ray trajectory consists of segments inside the homogeneous cheese and segments 
inside voids which are glued together in the global coordinates ( |3,4| ). The final conditions 
of a cheese-segment determine the initial conditions for integration inside the following 
void-segment. 

3.3 Void-segments 

The geodesic equations in the plane of the ray (eft = const) are 



d 2 t R'R' fdr\ 2 ■ fd9\ 2 , . 



d 2 r R' dt dr fR" E' \ ( dr\ 2 . _ R ( d6\ 2 , , 

d 2 9 RdtdO R' dr dO _ d ( R2 d9 \ _ 
d>? +2 RdXdX +2 ^dXdX ~ dX V d\)~ 0, (3 ' ' 

where A is the affine parameter along the geodesic. The null condition is the first integral 



1 + 2E \d\J \dXJ \dXJ 

To decrease the computational time, the total differential order of the system can be reduced 
by using first integrals: 

f = 4<< I + 1)» -(*/«)') -jj* (3-9) 



dX ^ R> v v 
dd L 

dX ~ ' 

where the last equation is an integration of ( |3.7|) and L is the integration constant analogous 
to angular momentum (naturally conserved in spherical symmetry) . The first two equations 



above are a rewrite of equation (3.5) in which dr/dX was expressed from the null condition. 
The third equation is a square root of the null condition. Its integration has to be stopped 
at the ray's turning point where dr/dX = and restarted after that. The sign in the third 
equation is determined by whether the ray approaches the turning point or goes away from 
it. 

One needs the initial conditions of integration at point F on Fig.||. There the ray exits 
the void forward in time or equivalently enters the void when traced back in time. The 
triangle ABF consists of physical lengths measured by the matter-comoving observer at 
point F for an infinitesimal ray segment. The physical length of the ray segment is dl = \dt\ 
because the locally measured speed of light is dl/\dt\ = 1, according to the fundamental 
axiom of General Relativity. The null condition (|3.8| ) shows that ABF is a right triangle 
(simply because the coordinates fl2.2| ) are orthogonal). The angle a is obtained using the 
end point of the previous cheese segment: 

a = &rccos(dir ■ FC) , (3.10) 

where dir shows the spatial direction of the ray traced back in time and the dot denotes 
scalar product in the global coordinates (jO]). It is easily calculated since the global 
coordinates of F and dir are the end conditions of the previous cheese-segment. One can 
write the equations 

dt = (z + l)dX (3.11) 

R' dr 

dt cosa (3-12) 
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Rd9 = dt sina . (3.13) 

All the functions must be evaluated at F where R = r v a{t) and E{r v ) = 0. After little 
algebra one gets: 

z(F) = z F (3.14) 

-&, + !) 
dr {zp + 1) cosa 

dX o,{tp) 
d6 (zp + 1) sina 



dX r v a{tp 
,d9_ 
dX 



L = —R -pr\F) = r v a(tF){zF + 1) sina, 



where ij? is the time at which the ray hits point F. The initial condition used for integration 
is the one for L. 



3.4 Cheese-segments 

The final point of the previous void segment, corresponding to the smallest A when tracing 
back in time, is denoted again with F. For a cheese-segment, the origin of the local 
spatial coordinates can be conveniently chosen to coincide with F since there is no center 
of symmetry analogous to the void center. Rays follow straight lines in the homogeneous 
cheese. The M-axis of the coordinates can be oriented along the ray, u = dir, which sets 
9 = and L = 0. The geodesic equations ( |3.9| ) simplify to 

l = -| (2 + 1 » 2 < 3 - 15 > 

dt 

Tx= z + 1 

dx 2 + 1 

dX = ~ a(t) ' 

where a(t) = (t/to) 2 / 3 is the EdS scale factor. The ray trajectory in global coordinates is 

x = F + x(A) dir . (3.16) 

The spatial vector dir gives the direction along the ray traced back in time. The minus 
sign in the third equation above means that x(X) increases when A decreases (going back 
in time), corresponding to going away from the point F in the direction of dir as it should 
be. The initial conditions are the obvious 

z(F) = z F , t(F)=t F , x(F) = 0. (3.17) 

The corresponding analytic solution is: 

BEl + M(l + ^) (3.18) 
dtp 

t(A) = t F n 3 / 5 

z(A) = (z F + l)u- 2/5 - 1 

X (A) = ^£(i_ u i/5). 
ap 

3.5 Tracing luminosity distance 

The area distance by definition is equal to \J A/£l s , where VL S is an infinitesimal solid angle 
subtended by a conical ray at the point source and A is the ray cross-sectional area at 
some distance from the source. If the ray has an opening angle f3 <C 1 at the source, then 
S7 S = 2n(l — cos/3) ~ 7r/3 2 . For rays with a negligible shear, which is the case considered 
in the present paper, the cross section at some distance is a circle of radius I. The angular 
diameter distance at that position is by definition 



equal to the area distance. The reciprocity theorem, first derived in 1933 by Etherington 
and popularized by Ellis [29], connects vq (original notation) - the area distance for a 
past directed ray emitted from the Earth observer towards the supernova, to tq (original 
notation) - the area distance for a future directed ray emitted from the supernova towards 
the Earth observer: 

r G = r (l + z) , (3.20) 

where z is the redshift between the supernova and the Earth. The luminosity distance to 
the supernova is [29] 

d L =r G (l + z) = r (l + z) 2 . (3.21) 

It is quite surprising that these relations are true in arbitrary spacetimes not only in 
homogeneous ones. 

The angular diameter distance along a past directed ray obeys, Jl( 

1 d 2 d A 1 <l'\'~\ 1 ,„,„ 2 



-R^Wk" - a 2 , (3.22) 



R^W = 8irG(T fU/ k^k v --k-k) = 8TrGT 00 (k ) 2 . (3.23) 



d A d\ 2 y/A d\ 2 2 

where a is the ray shear, R^ v is the Ricci tensor and k^ = dx^/dX is the null tangent 
vector of the ray. The inverted Einstein equation R^ v = 8t:G{T^ u — g^ u T/2), where T^ v 
is the energy-momentum tensor and T is its trace, can be used to calculate the righthand 
side of the above equation: 

T. 
2 

The last equality follows from k ■ k = and the fact Too = p is the only non-zero component 
for a presureless dust in the comoving coordinates ( |2.2[) . Finally, k° = dt / dX = z + 1 and 
the angular diameter equation is 

h^ = ~^ Gp{x){z{x) +l)2 ~ a2 - (3 - 24) 

The initial conditions at the Earth observer are 

dA (0) = 0, ^)=-l. (3.25) 
dX 

The last equality is a consequence of the relation dd A ~ dl, where the physical distance 
along the ray close to Earth is dl = \dt\ = \dX\ due to the local speed of light c = \dl/dt\ = 1 
and the normalization choice (3.2). The negative sign reflects the choice that A decreases 
tracing the ray back in time, while the distance d A increases. 

Equation ( 3.24 ) has to be integrated along with the geodesic equations on each segment 
of the trajectory. The density p(X) is calculated from the first equation in ( |2.4| ) in which 
r = r(X) and t = t(X) are the ones obtained from integrating the geodesic equations. Once 
the angular diameter distance d A = is known, the luminosity distance is obtained from 
equation (3.21). 



The shear is generated by inhomogeneous distribution of matter, [10]: 



where L is the angular momentum integration constant in (|3.9| ) and 9 e = (l/y/A)dvA/d\ = 
(l/djCjddA/dX is the expansion of the ray bundle. The righthand side is zero in the 
homogeneous cheese and the solution is decaying, a = const d^ 2 . In an inhomogeneous 
void, one can estimate L 2 /R 2 = (d6/dX) 2 R 2 ~ (ir/2r v ) 2 r 2 ~ 1, the bracketed term in 



(3.26) is ~ 1 within the void wall and zero outside and the solution is a ~ AirGpAX with 
the wall thickness A A ~ 5 Mpc. The ratio of the two terms in equation ( 3.24Q can be 
estimated as a 2 / (AtrG p) ~ AitG p AX 2 . The current homogeneous density of the model 



p(to) ~ 10 8 for which AttG p AX 2 ~ 10 6 and the shear can be safely neglected in ( 3.24| ). 



It plays a significant role only close to mass concentrations like black holes, [10|. 



4. Distance modulus averaging for a single void 

Discussing a universe with a single void lays the foundation for understanding the case 
of many voids. The Earth observer will always be at the origin of the global spatial 
coordinates. The void center is placed on the X cLXIS, clS shown on FigJq. The past- 
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Figure 6: Single void on the X axis pene- 
trated by past-directed rays emanating from 
the Earth observer at the origin O. The min- 
imal cone for averaging is defined by opening 
angle /3 max . 
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Figure 7: Solid line - redshift z(X) for a 
ray emitted at 2 angle with the X axis and 
penetrating a void centered at X = 60 Mpc. 
Dashed line - redshift z(X) for the homoge- 
neous EdS. 



directed light rays look unremarkably straight, at least around time to. This is due to 
the chosen spatial coordinates inside the void ( |2.7D , which are approximately Newtonian 
around time to- The density inhomogeneities are not sufficiently large to cause a noticeable 
light deflection in Newtonian coordinates. 

4.1 Redshift surfaces 

The computed redshift along a given ray is shown on Fig.^. The two "bumps" on 
the plot are where the ray encounters the front and the back of the void. The peculiar 
velocity in the front is towards the Earth observer and the redshift is lower; the velocity 
in the back is pointed away and the redshift is higher. The matter in the EdS cheese has 
a zero peculiar velocity - it does not feel the presence of the void thanks to the matching 



conditions (2.6). As a consequence, the redshift outside the void is indistinguishable from 
the one of homogeneous EdS as Fig.0 indicates. The inhomogeneities, through which the 
ray passes, do cause a minute change in the redshift but it is a second order relativistic 




effect [15|. Primary changes in the redshift result from peculiar velocities of the supernova 
sources inside the void. 

Redshifts in the bumps on Fig]?] cor- 
respond to three A positions along the ray. 
For a universe filled with many voids, some 
rays encounter a double bump when they 
exit a void wall and enter another void 
nearby. Redshifts in such double-bumps 
can correspond to five A positions. 

All supernovas with the same redshift 
lie on a spatial redshift surface. A few 
of those are shown on Fig^j - the actual 
surfaces are obtained by rotation around 
the symmetry X axis. The numbers that 
label the surfaces are the corresponding 
redshifts multiplied by 10 4 . Surfaces 70, 
75, 203 and 212 consist of three sheets for 
certain directions - three points giving the 
same redshift along a ray, corresponding 
to a "bump" on Fig.[7[ The evolution of 
the triple surface is traced on Fig|| as 
the redshift increases to the right, two of 

the sheets of surface 70 merge on the X axis and open up into surface 75 which gradually 
straightens out into surface 128 somewhat before the middle of the void. The reverse 
metamorphosis is seen continuing to the right. Surface 128 folds back into surface 203 
which pinches off into the lens-like and open sheets of surface 212. At some distance before 
and after the void, the redshift surfaces 58 and 235 are just trivial spheres centered on the 
Earth observer. The surfaces become more symmetric with respect to the middle of the 
void as its distance to the origin increases. 

4.2 Averaging by solid angle and by mass 

The supernova data displays not the luminosity distance di itself but a logarithmic measure 
of it, the distance modulus: fi = 5 log 10 (c2r,/10f?c) . The goal of this paper is to compare 
the distance moduli of the inhomogeneous universe and the homogeneous background EdS: 



Figure 8: Redshift surfaces for a void centered at 
60 Mpc. The bold numbers show z x 10 4 . Some 
surfaces have multiple sheets. The void boundary 
is shown by a dashed line. 
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where di{z) is the inhomogeneous luminosity distance (numerically computed by ray trac- 
ing) at redshift z and df / dS (z) = 3to(l + z — \f\ + z) is the homogeneous one. For small 
corrections AdL = dh~ df dS <C df dS , the logarithm can be expanded: 



(4.2) 



showing that to first order the change in the distance modulus is given by the fractional 
change in the luminosity distance. 

In an inhomogeneous universe, di(z) depends also on the ray direction and so does 
Afj,(z) . The correction, averaged over all directions, is obtained by a weighted integration: 

{A ^ )) = fMp^)^l t (4 .3) 
J dW(dir) 

where dW{dir) is the probabilistic weight assigned to the infinitesimal element of the 
spatial surface or volume over which the averaging is performed. The element is seen in 
the direction dir from Earth. The probability for observing a supernova emission from 
that element is dP = dW(dir) / J dW(dir). 

The previous works assume the probability is proportional to the observational solid 
angle on Earth, dW = dQ, and the averaging is over a redshift surface. That is reasonable 
if all light sources are in the cheese where the density is homogeneous and the probability 
for a supernova emission is uniform. For practical calculations, the full Air solid angle of 
the Earth observer is broken down into small solid angles AS7j. A ray i with direction 
inside AO, is chosen as its representative. The averaging integral is approximated by the 
sum 

V • An,- AQa 

where A/ij is the distance modulus correction (with respect to the homogeneous EdS) 
evaluated at the points where ray i pierces the sheets of the redshift surface z. One ray 
will usually have several such points and a corresponding number of repetitions of AOj. 
As a result, the total solid angle will be Y^i A£l{ > An if the sheets are more than one 
anywhere, which often happens in an inhomogeneous universe. The averaging by solid 
angle becomes nonphysical if the redshift surface cuts through a void since a significant 
probability of supernova emission will be assigned to the void interior which is almost 
empty. The results of this type of averaging will be presented only for reference purposes 
to show that the nonzero corrections are not artifacts of the more physically appropriate 
averaging described below. 

The present paper assumes that the probability for observing a supernova emission 
from a given comoving volume is proportional to the rest mass inside it, dW = dm, in 
agreement with the astronomical views on supernova bias. The rest mass inside the volume 
specified by an observational solid angle on Earth d£l, and sandwiched between redshift 
surfaces z and z + dz is 

dm = Y pdAdl = Y P<Pa ^ dQdz, (4.5) 
^ ^ \dz/d\\ 



where the redshift surface area cut off by rays with directions within dQ is dA = d\ dQ 
(from (|3.19|) ) and the physical distance between the surfaces is dl = \dt\ = (z + 1) \d\\ = 
(z + 1) \dz\/ \dz/d\\. The sum is over all sheets of the redshift surface, if many. The 



averaging integral is discretized analogously to the previous averaging procedure: 
E* Arm _ Ei A ^ ^ 4, (* + l)AzA^ / |cte/dA|i 



<AM*)> 
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ZiPi d% (z + l)AzAtti I \dz/d\\i 



(4.6) 



The sum runs over all rays and all points where ray % pierces the sheets of redshift surface 
z. 

4.3 Maximal average correction in the minimal cone: numerical results 

The real universe does not contain vast homogeneous regions in between the voids. Such 
regions have = and will damp the correction coming from the inhomogeneous void 
when averaging is performed. The largest possible average (A[i(z)) is obtained when the 
amount of cheese is minimal. With that in mind, the averaging for a single void is restricted 
to rays inside the minimal cone with an opening angle f3 ma x on Fig|| The minimal cone 
contains the void and a minimal amount of cheese. The obtained average corresponds 
to the best case scenario when all the voids in that range are exactly at the same radial 
distance from the Earth observer, forming a spherical shell, and their average corrections 
add up constructively to the total average. 

The void is centered at distance d = 60 Mpc on the X axis. Due to the axial symmetry, 
it is sufficient to shoot rays only in the upper half of the XY plane and the full picture is 
obtained by rotation around the X axis. The angle from zero to f3 max = &ics'm(r v /d) = 30 ° 
was divided amongst 3000 rays separated by A/3 = 0.01°. Ray i represents the solid angle 
Af2j = 27r(cos(/3j — A/3) — cos(/3j)) between the cones with opening angles /3j — A/3 and 
/3j. The redshifts for averaging were 500 values equally spaced on the interval 0.001 < 
z < 0.025. For each redshift, the points where the rays intersect the redshift surface were 
found and the average correction (A/j,(z)) was calculated according to ( [4.4D and (|4.6| ). The 
results are shown on Fig.||. 

Although the solid angle average is not 
physically sensible for sources inside the 
void, it has the same general behavior as 
the mass average. This shows the results 
of averaging are qualitatively independent 
of the chosen averaging procedure. The 
redshifts with positive average corrections 
correspond to surfaces before 128 on Fig.|8[ 
which are shifted forward with respect to 
the EdS surfaces, thus having a bigger d^. 
Surface 70 (redshift 0.0070) is the one with 
the maximal mass-averaged correction on 
Fig.^. Surfaces after 128 have negative 
corrections since they are shifted backwards 

with respect to the EdS redshift spheres. The smaller magnitude of the negative correc- 
tions compared to the positive ones is caused by the bigger amount of cheese in the minimal 
cone at higher redshifts which damps the average towards zero. The positive and negative 
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Figure 9: Correction (Afi(z)) averaged by solid 
angle (dashed) and by mass (solid) in the minimal 
cone of a void centered at 60 Mpc. 
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corrections become more symmetric with each other for voids at larger distances from the 
origin - they are completely symmetric in the limit of an infinite distance. 

The small spikes appearing on Fig.0 and later figures are errors from the numerical 



discretization (4.6) of the integrals in ( f4.3|) . Notice they appear at redshifts in the proximity 
of surfaces z = 0.0084 and z = 0.0186 on Fig.|8] which have a segment pointing straight to 
the observer. It is hard to integrate over such a segment with finite number of rays in (|4.6j). 
Increasing the number of rays decreases the oscillations but becomes time consuming. The 
reader should simply ignore the oscillations. 
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Figure 10: Total solid angle for each red- 
shift surface used for averaging. 



Figure 11: Total rest mass for each redshift 
surface used for averaging. 



The observer may want to bin the data and find the average correction, [A/i], over 
a redshift interval. Naively, one can sum up all the 500 corrections on Fig|] and divide 
by their number to obtain [A/i] no j Ve = 0.014 for the angle- averaged and [Afj] na i ve = 0.010 
for the mass averaged correction over the whole plotted interval. More appropriately, one 
has to weigh each (Afj,(z)) by the total weight of its redshift surface z, shown on Fig.|l0| 
and Fig. 11. For a redshift surface with a single sheet, Qtat = 2tt(1 — cos f3 max ) = 0.84 srad, 



which is the horizontal line on Fig. 10. The two peaks on that figure correspond to the front 



and back of the void where redshift surfaces have several sheets leading to overcounting 



Afij. The weights on Fig.|ll| are increasing because there is more mass in the minimal 
cone at larger redshift. The two peaks in the front and back of the void are again due to 
overcounting caused by multi-sheet redshift surfaces there. The weighted averages over the 
redshift interval shown on Fig. ^ are [A/i] = 0.016 for angle averaging and [A/i] = —0.006 for 
mass averaging. The interval averages [A/i] are significantly smaller in magnitude than the 
single-redshift averages (A/i(z)) on Fig.^. That results from the tendency of the positive 
and negative corrections (A/i(z)) from the front and back of the void to cancel out. The 
cancellation is imperfect close to the Earth observer since the redshift surfaces are not 
exactly symmetric with respect to the middle of the void. For voids at larger distances, the 
symmetry improves and so is the cancellation. A perfect symmetry is achieved when the 
void is at infinite distance from the origin and the rays run parallel to the X axis. Then 
the interval average would be exactly zero. If one wants to avoid the corrections shown on 
Fig.0, the data must be averaged over intervals bigger than a single void. 



The interval averages calculated above implicitly assume that the supernova detection 
efficiency is constant throughout the redshift interval. That is usually not true in obser- 
vational astronomy since supernovas at higher redshift are harder to detect and less in 
number. The detection efficiency would be especially non constant for big redshift inter- 
vals/big voids. That would bias the interval average [Afj] towards positive values. 



4.4 Maximal average kinematic correction vs. photon number conservation 



The argument [|T^] based on the gravitational lensing conservation of the total photon 
flux states that the luminosity distance, averaged over all directions in an inhomogeneous 
universe, is the same as in the homogeneous counterpart. The inhomogeneities induce 
focusing/defocusing of light but the argument ascertains these are random along different 
lines of sight and cancel out in the average. In reality, they are correlated close to the 
observer and the average does not have to be zero. 

The argument in [12] is implicitly staged in the Newtonian coordinates (|2.8| ). It assumes 
that the redshift surface in the inhomogeneous universe remains spherical and at the same o 
coordinate distance from Earth as in the homogeneous counterpart. The analysis in the 
previous sections shows this is not the case because the matter peculiar motion can shift 
and bend the redshift surface as shown on Fig|| If the redshift surface shifts along a 
ray by a comoving distance As in Newtonian coordinates, the kinematic correction to the 
luminosity distance is 

Ad L ^(l + z)As, (4.7) 



Equation (4^) shows the distance modulus correction Afj, is proportional to the fractional 
change in the luminosity distance Ad^/di ~ As/s (s is the comoving distance to the 
redshift surface along the ray). It will decrease as A[i ~ 1/s, because the kinematic shift 
As is bounded from above by the void comoving radius. For distances much bigger than 
a void radius, the photon conservation argument will apply but only in the approximate 
sense that with the distance s increasing, the redshift surface is getting closer to spherical 
and the fractional corrections are getting smaller. 

From the linear perturbation formula ( |2.14| ), the shift of the maximal correction surface 
70 on FigJ|| is As max oc v max /H Ft ar^\5\/2> = r_|<5|/3(l + z), where r_ = 25Mpc is the 
radius of the underdense void interior. For models with a linear growing perturbation in 
an EdS background \5\ = |e>o| a, where \5o\ = 0.9 is the underdensity in the void interior at 
present time. The change in the luminosity distance is Adi ~ (1 + z)As max oc r_|<5o|/3(l + 
z). The corresponding maximal average correction for a given redshift is obtained from 
P): 

r?|(5o|r_ 
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There are additional corrections to di due to weak lensing or non-linear effects like Rees- 
Sciama. Neglecting them leads to a very good fit to the numerically calculated maximal 
average correction at low redshifts which proves the kinematic correction Q4.7| ) dominates o 
there. The other effects become apparent at high redshifts where the kinematic correction 
has already decayed and for big voids in a deep nonlinear regime or with a decaying mode. 
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Figure 12: Solid curves - mass-averaged (/j, — /i EdS } in the minimal cone of a void centered at 
X =(60, 120, 180, 240, 300, 360, 420, 480, 540, 600) Mpc. Dashed curves - maximal average 



correction (4.8) for j] = 0.20 (upper) and T) = —0.18 (lower). Dotted curve shows n ACDM ~ yp dS . 



where df dS (x 



\z) = 3t (l + z- VT + z) is the homogeneous luminosity distance. The pa- 
rameter rj represents the ratio of the maximal average correction to the naive maximal 
correction corresponding to the maximal peculiar velocity. The value of n is not derivable 
from linear perturbation theory but requires a numerical calculation. It depends on the 
details of the model (peculiar velocities) and the specifics of the chosen averaging proce- 
dure (mass distribution). It was found numerically that rj does not depend on redshift for 
the models with a pure growing mode considered here. Models that contain a mixture of 
growing and decaying mode have a more complicated dynamics; the redshift surface shapes 
have a strong time dependence and rj ^ const. 

For small corrections, the logarithm in (£Oj) can be expanded, showing the correction 
decays like r_/s or in other words oc 1/N V where N v = s/(2r v ) is the number of void 
diameters that equal distance s. 



Fig. 12 is a composite plot showing the distance modulus correction averaged by mass 



in the minimal cone of a void placed at various distances from the Earth observer. The void 
positions do not intersect each other but the corrections overlap since a void can influence 
a nearby redshift surface that is outside - surfaces 70 and 212 on Fig.|| gained extra sheets 
because of the void. The positive and negative corrections get more symmetric at bigger 



distances from the observer. The dashed curves on Fig. 12 show the maximal average cor- 
rection estimated by (|4,8|). The upper curve, enveloping the maximal positive corrections, 
corresponds to w = +0.20 and the lower curve, describing the negative corrections, is for 
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rj = —0.18. At small redshifts, the voids enter nonlinear regime and the correction is higher 
than the linear estimate ( [4.8[ ). The dotted curve on Fig. 12 corresponds to the standard 
ACDM with £l m = 0.3 and Q\ = 0.7 - obviously the average correction cannot substitute 
for dark energy at redshifts larger than a void diameter from the observer. 

The simple estimate (|4.8| ) of the ideal scenario maximal correction does not take into 
account the presence of an excessive amount of cheese in-between the voids nor the fact 
that the corrections from randomized voids often add destructively as will be seen later. 
The first effect is an artifact of the Swiss-cheese model - the real universe does not contain 
vast homogeneous areas between the voids. Both effects will damp the averaged corrections 



of the model and they will go to zero with redshift faster than shown on Fig. 12 



5. Cumulative distance modulus correction by aligned voids 



The cumulative distance modulus correc- 
tion for a ray at 0.02 degrees angle piercing a 
string of 60 voids (r v = 30) aligned along the 
X axis is shown on Fig.|l3[ The voids touch 
but do not intersect. Each dot shows the cor- 
rection right after the ray exits a void. The 
correction inside the voids is not shown - it is 
multi-valued for some redshifts due to multi- 
sheet redshift surfaces. Including the ray shear 



in the distance tracing equation (|3.24|) does not 
influence the result. A correction of a similar 
magnitude at large redshifts was obtained in 
Q but for much larger and significantly non- 
linear (density contrast in the void wall 5 (to) 
effect was ascribed to the cumulative gravitational lensing from the underdense void in- 
teriors pi and shown to go to zero when averaged over randomized impact parameters. 



Figure 13: Correction A// for a ray at 0.02 
degrees with the X axis penetrating a string of 
60 voids centered at X = 2r v i with i = 1...60. 

= 28 0) voids of radius 350 Mpc. The 



Fig. 13 demonstrates that a large correction is possible for small voids too. The cumulative 



correction at a given redshift does not depend strongly on the size or the nonlinearity of 
the voids but only on how empty they are along the ray. 

One of the main questions addressed in this paper is whether the significant cumulative 
correction at big redshifts persists in the distance modulus when it is averaged over all 
directions in a lattice of voids. The answer is no, as shown in the next section, and is 
obtained without any approximations or assumptions unlike the previous studies |11]. The 
more surprising fact is that randomization of the voids is not necessary - even regular 
void lattices like the simple cubic have a vanishing average correction at high redshifts. 
Of course regular lattices will still show large cumulative corrections in certain directions, 



as [30 1 demonstrates, but the required void alignment is extremely improbable in the real 
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6. Averaging in lattices of r v = 30 Mpc voids 



The real universe is permeated by a lattice of voids. Its effect on the averaged distance 
modulus is calculated by shooting past-directed rays from the Earth observer in directions 
covering the full solid angle. The luminosity distances at selected redshifts are averaged 
over all directions via fl4.6j). 

A lattice of 1691 voids was obtained by randomly generating void center coordinates 
(x, y, z) on the interval —390 < x,y, z < 390 Mpc. The voids can touch but not intersect 
each other. They are not allowed to contain the Earth observer or get closer than 1.5 r v to 
her since that would lead to too big (A//) corrections. The closest void center is at distance 
d = 46.9 Mpc from Earth. The amount of unwanted cheese between the voids depends 
on the packing efficiency of the lattice, defined as the total volume fraction that is inside 
voids. The maximal packing efficiency of identical spheres is ir/y/l8 ~ 0.74 achieved in the 
face-centered cubic and the hexagonal close-packed lattices. The generated random lattice 
has a packing efficiency of 0.35 implying too much cheese. Unfortunately, higher packing 
cannot be achieved by the simple random process described here. 

For comparison, the distance modulus correction was also averaged in a simple cubic 
void lattice with a cube cell size of 2r v and a packing efficiency of 7r/6 ~ 0.52. In this case, 
one can average only over the first octant, the others being identical, thus reducing the 
computational work 8 times. The closest void in that lattice is at distance d = v3r„ « 51.96 
Mpc from the Earth observer. 

For the random void lattice, the total observational solid angle of An was divided 
approximately equally between 208024 directions/rays. The resolution in the observational 
polar angles was d6 = d<ft = 2/257 i.e. the rays could resolve a detail of a comoving size 2 
Mpc at a comoving distance 257 Mpc (corresponding to redshift z = 0.063). The luminosity 
distance was traced along each ray and its values corresponding to 200 redshift points 
equally spaced on the interval < z < 0.08 were extracted. Those values were used in (fOj) 
to calculate the average for each redshift point. The resulting mass-averaged correction is 
shown on Fig. 14. For the simple cubic lattice, the 4-7r/8 solid angle of the first octant was 
divided between 83599 directions/rays. The chosen resolution was dO = d<f> = 1.5/345 i.e. 
the rays could resolve a detail of size 1.5 Mpc at comoving distance 345 Mpc (corresponding 
to redshift z = 0.086). The redshift points were the same as for the random lattice and 
Fig. 15 shows the obtained average correction. The thin solid line shows the standard 
deviation of the corrections at each redshift <r[A/x](z) = y var[Afi] with the variance 
var [A/x] = (A// 2 ) - (A/z) 2 . The averages (...) over directions were calculated with the 
mass weights dW = dm. 

The average corrections on Fig. 14 and Fig. 15 are damped below the maximal average 
correction (4.8) given by the dashed curve. The first reason for that is that the two lattices 
have packing efficiencies lower than the one in the minimal cone, which approaches ~ 2/3 at 
distances bigger than several radii. The bigger amount of cheese with zero peculiar velocity 
in the lattices compared to the minimal cone damps the average correction towards zero. 
The much more serious reason for the damping is that whenever a void front and another 
void back are at the same distance from the Earth, their corrections will tend to cancel 
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Figure 14: Thick solid curve - mass- 
averaged (A/x) for the random void lattice. 
Dashed curve - maximal average correction 
for r\ — 0.20. Thin solid curve - standard 
deviation of A/z at each redshift. 
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Figure 16: Distribution of voids in radial 
distance (Mpc) for the random lattice. 
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Figure 15: Same as the previous figure but 
for the simple cubic lattice. 
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Figure 17: Same as the previous figure but 
for the simple cubic lattice. 



out in the total average. The astonishing efficiency of this cancellation damping, especially 
in the random lattice, at redshifts even as small as a void diameter is puzzling. This can 
be understood visually by looking at Fig. 16 and Fig.[l7] which show the void distribution 
in radial distance for the random and simple cubic lattices. Obviously, the front-back 
cancellation is more severe in the random lattice and at higher redshifts since there is 
more volume available there and correspondingly more voids per redshift interval. Quite 
surprisingly, Fig. 15 demonstrates that the cancellation occurs very vigorously also in a 
regular void lattice and even there the cumulative corrections from the previous section do 
not survive in the directional average. 

The voids in the real universe are not exactly spherical which allows for a tighter 
packing and more correlated radial void positions close to the observer than in the latices 
considered here. The expected damping would be between that of the random and the cubic 
lattice which in this way serve as lower and upper estimators of the average correction. 
Further discussion of the numerical results is given in the summary and conclusions section. 



7. Averaging in a lattice of r v = 300 Mpc voids 

A natural question to ask is whether bigger voids would produce a proportionally bigger 
average correction than the one obtained so far. 

The current astronomical data suggests that there are structures of size significantly 
larger than 30 Mpc. Galaxy count studies imply the Southern local void may have a 
radius of 200 Mpc and a density contrast of —0.25. The size and underdensity of the largest 



voids is restricted by the imprint they would leave on the CMB passing through them. A 



recent paper [32] found such an imprint by averaging the CMB temperature around voids 
at a mean redshift z v = 0.53. Initially, the signal was attributed to the late-time integrated 
Sachs- Wolfe effect. Since a significant linear Sachs- Wolfe effect is absent in a flat matter- 
only universe (see references in |33[| ), the signal was presented as a direct manifestation of 
dark energy. However, two subsequent papers [|33|, [34| showed that the measured signal is 
orders of magnitude too large to be explained with the Sachs- Wolfe effect in ACDM. 

The matter density of the voids in this section is chosen to reproduce roughly the 



average imprint on CMB, shown as a dashed line on Fig. 9 in [33]. The necessary matter 
density, shown on Fig. 18, is 



p(r,t ) = p(t ) 



A 1 + A 2 tanh[0.2(r - 200)] - A 3 tanh[0.3(r - 290)] 



, r < 300 
, r > 300 



(7.1) 



with the coefficients (A x , A 2 , A 3 ) = (0.764715, 0.349972, 0.115257), determined by the 
mass-compensation condition (|2.6j ) and density continuity. A more precise match to the 



curve in 1 32] would require a density function with more parameters and is not the goal 
of the present paper. The bang time chosen for the model is ts(r) = 0, implying a pure 
growing mode. The homogeneous "cheese" between the voids is the same Einstein-de Sitter 
metric that was used for the small voids. 



P (r, t ) I p(t a ) 
1.2; 
1.0; 



. 6 
. 4 
0.2 



50 100 150 200 250 300 350 



r, Mpc 



4 6 8 10 12 

angle (degrees) 



Figure 18: Density ratio p/p at time to for Figure 19: Change in the CMB tcmper- 
a void of radius r v — 300 Mpc. ature [p,K) produced by a void of radius 

r v = 300 Mpc centered at redshift z v = 0.53. 



Fig.19 shows the effect of such a void, centered at redshift z v = 0.53 (1633 Mpc in 
EdS), on the CMB temperature. It is calculated by shooting rays as shown on Figj6| 
and propagating them to some redshift zf and time tF outside the void. The redshift is 

3/2 

propagated further to the time of last scattering tis = to a ls us i n S the analytic solution 
d3.18[ ): zls = ( Z F + 1)(^fAls) 2//3 — 1- The fractional change in the CMB temperature 

1099 is the unperturbed redshift 



)/z° LS , where z° LS 



is AT/T = -Az/z = -(z LS ~ z° LS 
corresponding to clls = 1/1100. 

The effect shown on Fig.|l9| and in [33] is orders of magnitudes smaller than the one 
calculated in [30|. The large AT/T in [30] is due to: (1) a cumulative effect along special 



directions from aligned voids in a regular lattice, (2) a significant gravitational blueshift 
of the rays starting from inside voids at the initial time of the model when the density 
contrast in the void interior is practically 5 (Unit) = ~~ 1> (3) possibly effects produced 
by the significant non-linearity of the chosen void model || having a density contrast in 
the void wall 5 (to) = 28 pi]] . In the real universe, special alignment of more than 2 
- 3 voids to produce a cumulative effect on CMB has a vanishing probability. The CMB 
radiation starts its journey towards us when the dark matter density contrast is of the order 
\S\ ~ 10~ 3 . That produces a gravitational blueshift much smaller than the one exhibited 
in |3(J • These considerations invalidate the claim in |3(J that voids of radius larger than 
35 Mpc are excluded because they modify significantly the large-scale part of the CMB 
temperature spectrum. Such voids are actually found |22| in the SDSS data. 

The cumulative correction /i — fj, E along a ray at 0.02 degrees with the X axis, 
penetrating a string of 6 non-intersecting voids centered at X = r v (2i — 1), i = 1...6 is 
A/i = 0.11 at redshift z = 1.97. This is about twice as low as the correction for the small 
voids (A/j = 0.27), see Fig.|b|. The density contrast in the interior of the bigger voids is also 
about twice as low as for the small r v = 30 Mpc voids which indicates that the cumulative 
effect at a given redshift does not depend on the void size but on how empty the voids 
are along the ray. Analogously to the case of small voids, the cumulative correction is 
completely destroyed in the averaging over all directions, even for the regular simple cubic 
lattice. 




o 
o 

ro 
o 
o 
<x> 

o 
o 
o 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 



Figure 20: Solid curves - mass-averaged 
(/i — ix Eds ) in the minimal cone of a void of 
radius r v — 300 Mpc centered at X =(400, 
800, 1200, 1600, 2000, 2400) Mpc. Dashed 
curves - maximal average correction for r\ = 
0.22 (upper) and 77 = -0.20 (lower). Dotted 
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Figure 21: Thick solid curve - mass- 
averaged (A/x) for the simple cubic void lat- 
tice. Dashed curve - maximal average cor- 
rection for rj = 0.22 and r v — 300 Mpc. 
Thin solid curve - standard deviation of A/j 
at each redshift. 



Fig. 20 shows the correction averaged in the minimal cone for several positions of a 



single void. It is counterintuitive that the correction at distance r v is not bigger than the 
corresponding one for the small void. Its magnitude is determined by the ratio Adj^f v v oc 
v max/ r -- The explanation is that although the bigger voids have proportionally larger 
peculiar velocities, the linear formula ( |2.14| ) shows that the ratio v max /r^ is independent 



of the characteristic void size r_. The dashed curves on Fig. 20 are the maximal average 
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correction estimate fl4.8| ) for r_ = 200 Mpc and \5o\ = 0.47. The linear perturbation theory 
estimate (|4.8| ) envelops the maximal correction solid curves very well except at low redshifts 
where the voids become mildly nonlinear. It is a little surprising that approximately the 
same values of r/ from Fig.^ work on Fig.2C for an entirely different void model. That 



suggest some kind of universality of those values for the mass averaging procedure, at least 
for mass-compensated voids that have a pure growing perturbation mode, are at most 
mildly nonlinear, and have approximately constant densities in the void interior and walls. 



Fig. 21 shows the correction averaged over all directions in the simple cubic lattice. 
The random lattice was not used since it underestimates the average by having too much 
cancellation damping as was discussed in the case of small voids. To calculate the average 
(A/i), the simple cubic lattice used for the r v = 30 Mpc model was rescaled to the new void 
radius r v = 300 Mpc by multiplying the void center coordinates by 10. The 4 ir/8 solid angle 
of the first octant was divided between 51126 rays with a resolution dO = deft = 10/1800 
at distance 1800 Mpc (redshift 0.60). The redshift points for averaging were 200 and 
equidistant on the interval < z < 0.8. The average correction (thick solid line) is damped 
with respect to the maximal average correction fl4,8| ) (dashed line) due to cancellation 
between fronts and backs of different voids, analogously to the case of small voids. Further 
discussion of the numerical results continues in the summary and conclusions section. 

8. Summary and conclusions 

The paper studied how the distance modulus, averaged over all lines of sight in an inho- 
mogeneous universe, differs from its homogeneous counterpart. The inhomogeneities were 
represented by random and regular lattices of mass-compensated voids in Swiss-cheese 
models. Two void layouts were considered with a small (30 Mpc) and a big (300 Mpc) 
radius, the first observed in the SDSS data p3] and the second deduced from its imprint on 
CMB [ |32| |. The Earth observer was put in the cheese but the conclusions will not change 
significantly for an observer inside one of the voids. For the first time, the averaging of 
the distance modulus was widened to include the supernovas inside voids. That was made 
possible by assuming that the probability of supernova emission from a comoving volume 



is proportional to the rest mass in it. Unlike previous studies [12, 11], the average correc- 
tion to the distance modulus was calculated by exact numerical ray tracing without using 
assumptions about the ray impact parameters or perturbation theory approximations, and 
the result includes all linear and non-linear effects. 

The distance modulus correction, due to gravitational lensing, that accumulates along 
a ray crossing diametrically a sequence of aligned voids g [R], was calculated. The 



correction increases linearly and becomes quite significant at high redshifts, see Fig. 13. It is 
roughly proportional to the density contrast in the void interior but does not depend on the 
void radius. It was demonstrated that the distance modulus, averaged over all directions 
in a lattice of voids, does not show a cumulative correction at high redshifts. That is 



true not only in a random lattice (see Fig. 14) but also in a regular simple cubic one (see 



Fig. 15 and Fig. 21) implying that void randomization is not necessary for the destruction of 



the cumulative correction when averaging, contrary to popular belief. A large cumulative 



correction is still observed along special directions in a regular void lattice |30[], but the 
probability for the required void alignment is vanishing in the real universe. 

At low redshifts z, perturbation theory predicts that the line-of-sight peculiar ve- 
locities v r are capable of producing significant fluctuations in the luminosity distance 
AcIl(z) I di fa —v r /z [35, [36| (v r measured in speed of light units). That leads to a scatter 
in the local Hubble diagram AH /H = —Adi/di and, by ( f4.2[ ), to a non- vanishing average 
correction (Afi(z)) fa —2.17 {AH/ H) fa —2.17(v r {z)) / z, where the so called peculiar 
velocity monopole (v r (z)) is the line-of-sight peculiar velocity at redshift z, averaged over 
all directions. A measure of the fluctuation magnitude is its average within a sphere of 
radius R around the observer, [AH/H]r, which varies with location. Its cosmic variance 
over all possible observers/locations has been calculated previously in several cosmological 
models by N-body simulations |37| and by linear perturbation theory using the matter 
power spectrum |38], |39|. Such variances are employed to predict confidence intervals for 
the value of [AH / H\r measured by a random observer. Since the obtained cosmic variance 



of [AH/H]r decreases with R [38, p9|] , with the cosmic mean being zero, the likelihood to 
measure large average fluctuations decreases with redshift, as expected. 

The cosmic variance of the sphere-average [AH/H]r over all observers cannot be used 
to estimate the average over all directions (AH(z)/H) in the redshift bin z, experienced 
by a particular observer. The results in the present paper allow to calculate (AH(z)/H) 
by studying such an observer in a simple model of our local neighborhood. The averages 
and variances so obtained are not cosmic, thus are more applicable to our astronomical 
observations. They cannot be calculated from a given matter power spectrum but require 
a concrete numerical realization of an inhomogeneous universe. To the best knowledge 
of the author, this is the first time when the non-vanishing direction-averaged correction 
(Afi(z)) is calculated as a function of redshift for universes containing voids. The upper 



bound fljD|) for (A^(z)) was motivated using linear perturbation theory and was found 



to agree with the numerical calculations on Fig. 12 and Fig.^0[ The parameter r\ depends 
on the particular form of inhomogeneities and the choice of averaging procedure. The 
numerically obtained value r/ ~ 0.20 indicates that the maximal average correction is 20% 
of the naive correction corresponding to the maximal peculiar velocity. This r\ turned out 
to be a universal constant, approximately independent of the void size, for voids which 
have approximately constant densities in their interior and walls that are not in a deep 
nonlinear regime. Due to void randomization, it is expected that the average correction will 
be suppressed below the maximal 20% and will tend to zero. The efficiency of that process 
and at what redshift it sets in has not been studied before. The calculations revealed that 
the average correction within void lattices is indeed severely damped below the predicted 
maximal average correction due to cancelations between the fronts and backs of different 
voids. As figures 14, 15, and 21 demonstrate, the cancelation is surprisingly efficient at low 
redshifts even in regular lattices and the average correction drops below 0.01 mag after a 
single void diameter. Nevertheless, the average correction is not zero close to the observer, 
indicating that the implicit assumptions of the photon flux conservation argument stated 
in [12] do not apply at low redshifts. 

The nonzero (Afi(z)) is observable, provided there are enough supernovas N z in redshift 



bin z to reduce the sampling statistical error of the average: a[ (A/i(z)) ] = cr[Afi(z)]/ y/Nz, 
where a[Afi(z)] is shown as a thin solid curve on figures 14, [l5|, and ^l|. An underdense 
"Hubble bubble" around us with a Hubble parameter slightly higher than the global one 
(correspondingly (Afi(z)) < 0) has been argued for in [[!(], ^TJ. However, that claim was 
not confirmed in p2| , |43] ] and was shown to depend on the way the supernova magnitudes 
were extracted from the photometric data [44|. Local bubble or not, fluctuations in the 
distance modulus binned average begin to appear in the newest data sets, Fig. 20 in 



and Fig. 7 in g6|. 

Significant efforts in contemporary cosmology are aimed at measuring a possible time 
evolution in the dark energy equation of state. That would require fixing the Hubble 
constant at low redshifts to a 1% precision going and future low redshift 

surveys will boost the number of the observed nearby supernovas sufficiently to bring the 
sampling error of the average a[(AH/H)} « a[ (A/i(z)) ] / 2.17 below 1%. In contrast, 
the directional average itself (AH/H) « — (Afi(z)) /2.17, generated by the coherent 
peculiar motion, is not influenced by the supernova number and will degrade significantly 
the errors in measuring the dark energy time dependence [p5 



Consequently, the 



peculiar velocities need to be subtracted from the low redshift Hubble diagram. 

Several methods are utilized for that. The Local Group (LG) of galaxies moves as 



a whole with respect to CMB at a speed vlg, the so called peculiar velocity dipole |50 
51], The most common way to correct the Hubble diagram for that motion is to adjust 



the observed redshifts to an observer riding the LG barycenter. The luminosity distance 
fluctuations seen in that frame are Ad^/dL ~ —(v r — vlg)/ z [HH an d the velocity dipole 
will cancel out since the coherent velocity component of a nearby galaxy with respect to 
CMB is Vr° h ~ vlg- Not correcting the CMB redshifts to the LG frame is permissible for 
a large supernova sample that covers densely and uniformly all directions: the coherent 
dipole fluctuation seen in the CMB frame is Ad^/di ~ —v r /z — vlg cos(8) jz and 
that expression vanishes when averaged over the full solid angle. A further refinement 
is to correct the LG redshifts for infall velocities in the LG frame caused by the nearby 
superclusters represented by a crude linear multi-attractor model. That was done in the 
Hubble Key Project Q| which measured the Hubble constant to a 9% precision, but it 
was not very efficient at reducing the scatter in the low redshift Hubble diagram as seen 



on Fig. 4 of [52]. Accordingly, the philosophy of the Hubble Key Project was to extract 
the Hubble constant mainly from secondary distance indicators at high redshifts. A finer 
method to correct for peculiar velocities is to calculate them from the observed galaxy 
distribution utilizing linear perturbation theory and choosing a galaxy-dark matter biasing 
parameter that minimizes the scatter on the Hubble diagram. That technique does reduce 



the scatter [53] but using it to infer the Hubble parameter from low redshift data produces 
a controversially high value of Hq = 85 km Mpc~ l [|54|]. To avoid unreliable velocity 
corrections, many recent surveys [45, 47, ^] simply chose to include only objects above 
a certain redshift, z > z m i n , where the peculiar velocities are considered insignificant 
compared to the Hubble flow. The Hubble constant and the dark energy equation of 
state parameter w inferred from the low redshift data are sensitive to the choice of z, 
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the literature, ranging from 0.01 [55| up to 0.023 [47], which indicates that there is not a 
universally accepted prescription. 

Figures [l4|, [l5|, and 21 of the present paper allow to select z m i n readily. The required 
AH /H = 0.01 precision corresponds to A[i goa i = 0.0217. The sampling error of the average 
is a[A/i(z)}/\/N z , where a[A/j,(z)] is given by the thin solid curves on the plots and N z is 
the number of objects/supernovas in redshift bin z. Provided there is sufficient statistics, 
the sampling error will fall below A/j, goa i. In that case the limiting value z m in is determined 
by the average (A/j,(z)) itself - a natural choice is the redshift at which the maximal aver- 
age correction (|4.8[) (the dashed curves) drops below Afi goa i or the redshift corresponding 
to a void diameter (where {Afi(z)) ~ 0.01 due to front-back void cancellation), whichever 
is lower. In comparison, estimates based on the sphere-average [AH/H]r are more pes- 
simistic: its cosmic variance in the ACDM model drops below 1% at z m i n = 0.05 |3^] (linear 
estimate). If our neighborhood contains voids of the size observed on the CMB imprint 
]32|], the required z m i n would be much larger, as Fig. 21 indicates. 

o 
o 



Excluding the objects below z min reduces the size of the sample and proportionally 
increases the statistical error [48]. The present paper suggests it is possible to preserve 
the low redshift data, instead of throwing it away, by collapsing it into an interval average 
[A/j] over a redshift interval containing a void diameter. These averages turn out very 
close to zero: [A/j] = 0.003 (z = . . . 0.025 on Fig.|l4|), [An] = 0.007 (z = . . . 0.025 on 
Fig.15), [A^] = 0.002 (z = . . . 0.25 on Fig.^T|). That stems from the cancellation between 
positive and negative lobes in (Afi(z)) and the mass averaging giving higher weights to 
higher redshifts where {Afi(z)} is closer to zero. A vanishing [Afi] ~ means [fi] ~ [fi] EdS : 

Y^Wii^Zi)) n^Wi^iH^Zi), (8.1) 

i i 

where Wi is the mass weight assigned to redshift bin Zi, (fJ>(zi)} is the distance modulus 
in that bin averaged over all directions, and the sum is over all the bins in the redshift 
interval. The lefthand side of the formula is calculated from the observational data. The 
weight Wi is proportional to the rest mass inside the redshift bin which can be estimated 
as the corresponding mass in the homogeneous background model. Unlike a traditional 
Hubble diagram that weighs all redshift bins equally, here Wi oc zf Azi at low redshifts. 
The righthand side of ( |S.lD depends solely on the Hubble parameter Hq at low redshifts 
and on additional cosmological parameters of the background model at higher redshifts. It 
allows one to extract the value of Hq from the low redshift data. This method will become 
possible for future surveys that: (1) have a dense full-sky coverage to calculate the average 
correction over all directions (Afx(z)); and (2) contain a large number of objects in each 
redshift bin so that (A/j,(z)) is readily observable, not swamped by the sampling error. 

If the observer is inside a void, the average correction (Afi(z)) starts with a negative 
lobe unlike the so-far considered case of an outside observer. A simple illustration of that 
is an observer at the center of a void, which will measure An(z) = n{z) - n EdS {z) < o for 
supernovas inside the void since the matter there is expanding faster than the background 
and it takes a smaller distance di to achieve the same redshift. The interval average [A/i] 
over a void diameter would still be very small since the mass averaging scheme adopted 
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in this paper ascribes lower weights to low redshifts. Additional studies are needed to 
investigate the behavior of [A/i] in the presence of superclusters - the other type of major 
inhomogeneities encountered in the universe. A plausible guess is that [Afj] vanishes on a 
redshift interval encompassing a supercluster, justifying the validity of Q8,l| ) in that case 
as well. 
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